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ABSTRACT 

We report on multifrequency observations of the 7-ray emitting narrow-line Seyfert 
1 galaxy PKS 1502+036 performed from radio to 7-rays during 2008 August-2012 
November by Fermi-Large Area Telescope (LAT), Swift (X-Ray Telescope and Ul- 
traviolet/Optical Telescope), Owens Valley Radio Observatory, Very Long Baseline 
Array (VLBA) , and Very Large Array. No significant variability has been observed in 
7-rays, with 0.1-100 GeV flux that ranged between (3-7)xl0~^ ph cm~^ s""'^ using 
3-month time bins. The photon index of the LAT spectrum (F = 2.60 ± 0.06) and the 
apparent isotropic 7-ray luminosity (Lo.i_iooGeV = 7.8x10*^ erg s~^) over 51 months 
are typical of a flat spectrum radio quasar. The radio spectral variability and the 
one-sided structure, in addition to the observed 7-ray luminosity, suggest a relativistic 
jet with a high Dopplcr factor. In contrast to SBS 0846-f 513, the VLBA at 15 GHz 
did not observe supcrluminal motion for PKS 1502+036. Despite having the optical 
characteristics typical of a narrow-line Seyfert 1 galaxy, radio and 7-ray properties of 
PKS 1502+036 are found to be similar to those of a blazar at the low end of the black 
hole mass distribution for blazars. This is in agreement with what has been found in 
the case of the other 7-ray emitting narrow-line Seyfert 1 SBS 0846+513. 

Key words: galaxies: active - galaxies: nuclei - galaxies: Seyfert - galaxies: individ- 
ual: PKS 1502+036 - gamma-rays: general 



1 INTRODUCTION 

A small fraction of narrow-line Seyfert 1 (NLSyl) ga lax- 
ies are known to be radio loud l|Komossa et al.l 120061 ) . In 
these cases, the flat radio spectra and flux density vari- 
ability suRses t that several of them could hos t relativis- 
tic jets (e.g. IZhou et all l2003l : iDoi et all l2006t ). The de- 
tection by the Large Area Telescope (LAT) on-board the 
Fermi satellite of variable 7-ray emission from some radio- 
loud NLSyl galaxies revealed the presence of a possible 
third class of 7-ray emitting active galactic nuclei (AGN; 



lAbdo et al1l2009l : [P^mmando et al.l 1201231 ) in addition to 
blazars and radio galaxies, both ho sted in giant ellipti- 
cal galaxies (|Blandford fc Ree3 19781). NLSy ls are usually 
hosted in spiral galaxies (e.g. iDeo et al.ll2006l ) and the pres- 
ence of a relativistic jet in these sources seems to be in con- 
trast to the paradigm that the formati on of relativistic jets 
could happen in elliptical galaxies only (jBottcher fc DermeJ 
2002; Marschcr 2010). This discovery poses intriguing ques- 
tions on the nature of these objects, the formation and de- 
velopment of relativistic jets, the mechanisms of high-energy 
emission, the AGN Unification model and the evolution of 
radio-loud AGNs. 



* E-mail: filippo.dammando@fisica.unipg.it 



PKS 1502+036 has been classified as an NLSyl on ba- 
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sis of its optical characteristics: full width at half-maximum 
(FWHM) (H^) = 1082 ± 113 km s-\ [OIII]/H/3 ~ 1.1, 
and a strong Fe II bump. For this source, a radio loudness 
= 1549 was estimated l|Yuan et al.l 120081 '). R being de- 
fined as the ratio between the 1.4 GHz and 4400 A rest-frame 
flux densitiej3- This is one of the highest radio-loudness pa- 
rameters among the NLSyls. With its convex radio spec- 
trum peaking above 5 GHz, PKS 1502+036 was selected by 
iDallacasa et al.l (|2000l ) as part of the 'bright' high frequency 
peakers (HFPs) sample of young radio source candidates. 
However, its polarization properties, the flux density and 
spectral variability are typical of blazars and it was removed 
from the young radio source sample (|Orienti et al.l |2007| . 
|2008| ). Snapshot very long baseline interferometry (VLBI) 
observations in the 1990's showed a com pact source (on 
scale s of ~20 mas) with inverted spectrum (jDallacasa et al.l 
1 19981 ). 

In the 7-ray energy band PKS 1502+036 has been in- 
cluded in both the first and second Fermi-LAT source cat- 
alogu es (IFGL and 2FGL; lAbdo et all I2OI0I : iNolan et all 
l2012h . In the past, the source was not detected by the en- 
ergetic gamm a-ray experiment tele scope (EGRET) at E > 
100 MeV (see lHartman et al.|[T999l '). 

In this paper, we discuss the characteristics of 
PKS 1502+036, one of the 5 NLSyl detected by Fermi- 
LAT with high confidence, determined by means of new 
and archival radio-to-7-ray data. The paper is organized 
as follows. In Section 2, we report the LAT data analy- 
sis and results, while in Section 3, we present the result 
of the Swift observations performed during 2009-2012. Op- 
tical data collected by the Catalina Real-Time Transient 
Survey (CRTS) are reported in Section 4. Radio data col- 
lected by the Karl G. Jansky Very Large Array (VLA), the 
Very Long Baseline Array (VLBA) interferometers, and the 
40 m Owens Valley Radio Observatory (OVRO) single-dish 
telescope are presented and discussed in Section 5. In Sec- 
tion 6, we discuss the properties of the source and draw our 
conclusions. Throughout the paper, a A cold dark matter(A- 
CDM) cosmology with Ho = 71 km s"^ Mpc"\ = 0.73 
and fJm = 0.27 is adopted. The corresponding luminosity 
distance at z = 0.4088 (i.e. the source redshift) is d^ ~ 2217 
Mpc, and 1 arcsec corresponds to a projected distance of 
5.415 kpc. In the paper, the quoted uncertainties are given 
at the la level, unless otherwise stated, and the photon in- 
dices are parametrized as dN/dE oc with F = a+1 (a 
is the spectral index). 



2 FEi?M7-LAT DATA: SELECTION AND 
ANALYSIS 

The Fermi-hPs!T is a pair-conversion telescope operating 
from 20 MeV to > 300 GeV. It has a large peak effective 
area (~8000 cm^ for 1 GeV photons), an energy resolution 
of typically ~10%, and a field of view of about 2.4 sr with an 
angular resolution (68% containment angle) better than 1° 
for energies abov e 1 GeV. Further det ails about the Fermi- 
LAT are given in lAtwood et al.l (j2009l ). 

^ It is worth noting that this definitio n of rad i o-loud ness is re- 
lated to that presented in Kellerman n et al ] l| 19891 ). Re cm — 
f„(6cm)/f^(4400 A), by means of R1.4 GHz =1-9 Re cm- 
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Figure 1. Integrated flux light curve of PKS 1502+036 obtained 
in the 0.1-100 GeV energy range during 2008 August 4-2012 
November 4 with 3-month time bins. Arrows refer to 2a upper 
limits on the source flux. Upper limits are computed when TS < 
10. 

The LAT data reported in this paper were collected 
from 2008 August 4 (MJD 54682) to 2012 November 4 (MJD 
56235). During this time, the LAT instrument operated al- 
most entirely in survey mode. The analysis was performed 
with the ScienceTools software package version v9r27pl. 
The LAT data were extracted within a 10° region of inter- 
est centred at the location of PKS 1502+036. Only events 
belonging to the 'Source' class were used. In addition, a cut 
on the zenith angle (< 100°) was also applied to reduce 
contamination from the Earth limb 7-rays, which are pro- 
duced by cosmic rays interacting with the upper atmosphere. 
The spectral analysis was performed with the instrument re- 
sponse functions P7S0URCE_V6 using an unbinned maximum- 
likelihood method implemented in the Science tool gtlike. 
A Galactic diffuse emission model and isotropic component, 
which is the sum of an extragalactic and residual cosmic ray 
background were used to model the backgrouncQ. The nor- 
malizations of both components in the background model 
were allowed to vary freely during the spectral fitting. 

We evaluated the significance of the 7-ray signal from 
the sources by means of the maximum-likelihood test statis- 
tic TS — 2Alog(likelihood) between models with and with- 
out a point source at the position of PKS 1502+036 
iMattox et al.l ligod ). The source model used in gtlike 
includes all the point sources from the 2FGL catalogue 
that fall within 15° from PKS 1502+036. The spectra of 
these sources were parametrized by power-law (PL) func- 
tions, dN/dE oc {E/Eo)~'^, except for 2FGL J1504.3+1029 
and 2FGL J1512.8— 0906, for which we used a log-parabola 
(LP), dN/dE oc E/E-"-^ iog(E/Eo) jLandau et all Il986l : 
iMassaro et al.l '2004), as in the 2FGL catalogue. A first 
maximum-likelihood analysis was performed to remove from 

^ http://fermi.gsfc.nasa.gov/ssc/data/access/lat/Background 
Models.html 
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the model the sources having TS < 10 and/or the predicted 
number of counts based on the fitted model Npred < 3. A 
second maximum-likelihood analysis was performed on the 
updated source model. The fitting procedure has been per- 
formed with the sources within 10° from PKS 1502-1-036 
included with the normalization factors and the photon in- 
dices left as free parameters. For the sources located between 
10° and 15° from our target, we kept the normalization and 
the photon index fixed to the values of the 2FGL catalogue. 

The 7-ray point source localization by means of the 
gtfindsrc tool applied to the 7-rays extracted during the 
51 months of observation (MJD 54682-56235) results in RA 
= 226.27°, Dec. = 3.45° (J2000), with a 95% error circle 
radius of 0.07°, at an angular separation of 0.01° from the 
radio position of PKS 1502-1-036. This implies a strict spa- 
tial association with the radio coordinates of the NLSyl PKS 
1502+036. Integrating over the entire period considered, the 
fit yielded a TS — 305 in the 0.1-100 GeV energy range, with 
an integrated average flux of (4.0 ± 0.4) x 10~* photons cm^^ 
s~^ and a photon index of F = 2.60 ± 0.06. The correspond- 
ing apparent isotropic 7-ray luminosity is 7.8x10*^ erg s~^. 
In order to test for curvature in the 7-ray spectrum of PKS 
1502-1-036, an alternative spectral model to the PL, a LP 
was used for the fit. We obtain a spectral slope a = 2.45 ± 
0.11 at the reference energy Eq = 300 MeV, a curvature pa- 
rameter around the peak /3 = 0.10 ± 0.06, with a TS = 307 
and an integrated average flux of (3.7 ± 0.4) X 10"** photons 
cm~^ s~^. We used a likelihood ratio test to check the PL 
model (null hypothesis) against the LP model (alternative 
hypothesis). These values may be compared by defining the 
curvature test statistic TScurvc=TSLP-TSpL=2 correspond- 
ing to ~1.4(T, meaning that we have no statistical evidence 
of a curved spectral shape. 

Fig. [H shows the 7-ray light curve for the period 2008 
August 4-2012 November 4 using a PL model and 3-month 
time bins. For each time bin, the spectral parameters of PKS 
1502-1-036 and all sources within 10° from it were frozen to 
the value resulting from the likelihood analysis over the en- 
tire period. If TS < 10, the value of the fiuxes were replaced 
by the 2 a upper limits, cal culated using the profile method 
(see e.g. iNolan eraDl2012h . The systematic uncertainty in 
the fiux is energy dependent: it amounts to 10% at 100 MeV, 
decreasing to 5 % at 560 MeV, and incre asing again to 10% 
above 10 GeV (|Ackermann et al]|2012bl ). Despite integrat- 
ing over three months, the source was not always detected. 
We test for f lux variability calcu lating the variability index 
as defined in lNolan et"al] (|2012l '). We obtain a TSvar = 25.2 
for 16 degrees of freedom, below the value of 32 defined to 
identify variable sources at a 99% confidence level. Thus, 
no significant variability was detected in 7-rays, as reported 
also in the 2FGL catalogue over the first 24 mo nths of Fermi 
operations for this source (jNolan et alJl2012l ). On weekly 
time-scales, the source was detected with TS > 10 only a 
few times, with a peak flux of (18 ± 6) x 10~* photons cm~^ 
s"^ in the 0.1-100 GeV energy range during 2012 July 7-13. 

Leaving the photon index free to vary in the flts for 
three 3-month periods of highest activity (flux > 6xl0~* 
photons cm~^ s~^ in the 0.1-100 GeV energy range) and Eo 
flxed to 548 MeV, as in the 2FGL catalogue, the flts result 
in photon index values F = 2.57 ± 0.16, F = 2.44 ± 0.15 
and F = 2.54 ± 0.18, in chronological order. No conclusive 
evidence for spectral evolution has been detected. By means 



of the gtsrcprob tool, we estimated that the highest energy 
photon emitted by PKS 1502-1-036 (with probability > 80% 
of being associated with the source) was observed on 2009 
October 30 at a distance of 0.05° from PKS 1502-1-036 with 
an energy of 19.6 GeV. 



3 SWIFT DATA: ANALYSIS AND RESULTS 

The Swift satellite (|Gehrels et al.l |2004| ) performed six 
observations of PKS 1502+036 between 2009 July and 
2012 August. The observations were performed with all 
three onboard inst ruments: the X-ray Telescope (XRT; 
iBurrows et al.ll2005l, 0.2-10.0 keV), th e Ultraviolet/Optical 
Telescope rUVOT: lRo"ming et al.ll2005l. 170-600 nm) and the 
Burst Alert Telescope (BAT: iBarthelmv et al.ll2005l . 15-150 
keV). 

The hard X-ray flux of this source turned out to be be- 
low the sensitivity of the BAT instrument for such short ex- 
posures and therefore the data from this instrument will not 
be used. Moreover, the source was not present in the Swift 
BAT 70-month hard X-ray catalogue (|Baumgartner et al.l 
l2013h . 

The XRT data were processed with standard procedures 
(xrtpipeline vO.12.6), filtering, and screening criteria by 
using the HEAsof t package (v6.12). The data were collected 
in photon counting mode in all the observations. The source 
count rate was low (< 0.5 counts s~^); thus pile-up correc- 
tion was not required. The data collected during 2012 Au- 
gust 7 and 8 were summed in order to have enough statistics 
to obtain a good spectral fit. Source events were extracted 
from a circular region with a radius of 20 pixels (1 pixel 
~ 2.36"), while background events were extracted from a 
circular region with radius of 50 pixels far away from the 
source region. Ancillary response files were generated with 
xrtmkarf, and account for different extraction regions, vi- 
gnetting and point spread function corrections. We used the 
spectral redistribution matrices v013 in the Calibration data 
base maintained by HEASARC. Considering the low number 
of photons collected (< 200 counts) the spectra were re- 
binned with a minimum of 1 count per bin and we used 
Cash statistics (jCashl 0*9791 ). We have fitted the spectrum 
with an absorbed power-law using th e photoelectric absorp- 
tion model tbabs ( Wilm s et al.ll2000l ). with a neutral hydro- 
gen c olumn density fixed to its Galactic value (3. 93x10^" 
cm "^: lKalberla et al.ll2005l ). The results are reported in Ta- 
ble[l] In the past, the source has not been detected in X-rays 
during the ROSAT all-sky survey, with an upper limit of < 
2.1xl0~^^ erg cm ~^ s~^ in the 0.1-2.4 keV energy range 
(|Yuan et al.l 120081 ). PKS 1502+036 seems to have been ob- 
served by Swift/XBT in a relatively bright state on 2009 
July 25 in the 0.3-10 keV energy range. In addition a possi- 
ble moderate change of the photon index was observed dur- 
ing 2009-2012 together with an increase of the count rate 
leading up to 2012 August 7-8, but no conclusive evidence 
can be drawn on consideration of the large uncertainties. 

During the Swtft pointings, the UVOT instrument 
observed PKS 1502+036 in all its optical ( v, b and u) 
and UV (ml, m2 and w2) photometric bands (|Poole et al.l 
l2008l : iBreeveld et~al]|2010l ). We analysed the data using the 
uvotsource task included in the HEAsoft package. Source 
counts were extracted from a circular region of 5 arcsec ra- 
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Table 1. Log and fitting results of Siui/t/XRT observations of PKS 1502+036 using a power-law model with A'^h fixed to Galactic 
absorption. 



Date (UT) 


MJD 


Net exposure time 


Net count rate 


Photon index 


Flux 0.3-10 keV 






(sec) 


{xl0~3 cps) 


(r) 


(xlO"^"^ erg cm~^ s 


2009-July-25 


55037 


4681 


7.5 ± 1.3 


1.2 ±0.4 


4.9 ± 1.1 


2012-Apr-25 


56042 


4807 


7.7 ± 1.3 


1.7 ±0.4 


4.0 ±0.9 


2012-May-25 


56072 


4635 


8.0 ± 1.3 


1.9 ±0.4 


3.7 ± 0.9 


2012-June-25 


56103 


5142 


8.4 ± 1.3 


2.2 ±0.4 


3.5 ± 0.7 


2012-Aug-07/08 


56146/7 


4925 


10.8 ± 1.5 


2.2 ±0.3 


4.0 ±0.6 



''Observed flux 



dius centred on the source, while background counts were 
derived from a circular region of 10 arcsec radius in the 
source neighbourhood. We combined the images collected 
on 2012 August 7 and 8 in order to achieve better statistics. 
The observed magnitudes are reported in Table [2] Upper 
Umits are calculated using the UVOT photometric system 
when the analysis provided a significance of detection < 3cr. 

Magnitudes in optical filters are comparable within 
the uncertainties between the different Swift observa- 
tions, and compatible with the values observed in 
the past (^=1 8.62-18.82; _B=19.11-19.27; (7=18.58-18.68; 
I Wills fc LvM3I1978| ). Recently, optical intraday variability 
with amplitu de as large as ^10% has been reported for PKS 
1502+036 bv iPaliva et^ (|2013l ). No evident infrared in- 
traday variability has been found with Wtde-field Infrared 
Survey Explorer (WISE) data, but a variation of 0.1-0.2 
mag corrected for measurement errors during two epochs 
separa ted by about 180 d has been reported in I Jiang et al.l 
ll2012[) . 



4 CATALINA REAL-TIME TRANSIENT 
SURVEY 

The source has be en monitored by the CRTS 
||http://crts. caltech.edu iDrake et~lil ] |2009l : lDiorgovski et al.l 
I2OIII ). using the 0.68 m Schmidt telescope at Catalina Sta- 
tion, AZ, and an unfiltered CCD. The typical cadence is 
to obtain four exposures separated by 10 min in a given 
night; this may be repeated up to four times per lunation, 
over a period of ~6-7 months each year, while the field 
is observable. Photometry i s obtained using the sta ndard 
Source- Extractor package (jBertin fc Arnout3 Il996l ). and 
roughly calibrated to the V band in terms of the magnitude 
zero-point. During the CRTS monitoring, the source showed 
a variability amplitude of 1 mag, changing between 18.7 
and 17.7 mag. 



5 RADIO DATA: ANALYSIS AND RESULTS 
5.1 VLA and VLB A data 

Multi-epoch and multifrequency observations of 
PKS 1502+036 were carried out with both the VLA 
and VLBA between 1999 and 2007 as part of the monitor- 
ing c ampaign of the 'bright' HFP sample (jDallacasa et al.l 
I2OO0I ). Results of the observatio ns performed b etwee n 
1999 and 2003 were presented in iDallacasa et all (|20o3); 



iTinti et all (|2005l ): lOrienti et all (^2006l ■ l2007^ . To study the 
proper motion, we complemented our VLBA observations 
with the 15 GHz V LBA data from the MDJAVE programms0 
(jLister et al.ll2009l l and additional 15 GHz archival VLBA 
observations. Logs of VLA and VLBA observations are 
reported in Table [S] 

Simultaneous VLA observations were performed in L 
band (with the two intermediate frequencies, IFs, centred 
at 1.415 and 1.665 GHz), C band (with the two IFs cen- 
tred at 4.565 and 4.935 GHz), X band (with the two IFs 
centred at 8.085 and 8.465 GHz), U band (14.940 GHz), 
and K band (22.460 GHz). At each frequency the target 
source was observed for about 1 min, cycling through fre- 
quencies. During each run, the primary calibrator 3C 286 
was observed for about 3 min at each frequency. The data 
reduction was carried out following the standard procedures 
for the VLA implemented in the National Radio Astronomy 
Observatory (NRAO) Astronomical Image Processing Sys- 
tem (AIPS) package. The fiux density at each frequency was 
measured on the final image produced after a few phase-only 
self-calibration iterations. In the L band, it was generally 
necessary to image a few confusing sources falling within 
the primary beam. The target source is unresolved at any 
frequency even with the VLA A-configuration. During one 
epoch, strong radio frequency interference at 1.665 GHz was 
present, and a measurement of the flux density was not pos- 
sible. Uncertainties on the determination of the absolute flux 
density scale are dominated by amplitude errors, which are 
about 3% in L, C, and X bands, about 5% in U band, and 
about 10% in K band. 

Simultaneous VLBA observations were performed in 
2002 January at 15.3 and 22 GHz, and in 2006 July at 1.7 
GHz, 2.3 GHz (S band), 5.0 GHz, 8.4 GHz and 15.3 GHz. 
The correlation was performed at the VLBA hardware cor- 
relator in Socorro and the data reduction was carried out 
with the NRAO AIPS package. After the application of sys- 
tem temperatures and antenna gains, the amplitudes were 
checked using data on 4C 39.25 (J0927+3902), which was 
used also as bandpass calibrator. The error in the absolute 
flux density scale is generally within 3%-10%, being highest 
in value at the highest frequency. 

The final images were obtained after a number of self- 
calibration iterations. Amplitude self-calibration was ap- 
plied only once at the end of the process, using particular 



^ The MOJAVE data archive is maintained at 
|http://www. physics.purdue.edu/MO JAVE| 



Table 2. Results of the Swift/VVOT data for PKS 1502H 
detection < 3cr. 
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036. Upper limits are calculated when the analysis provided a significance of 



Date (UT) 


MJD 


V 


b 


u 


wl 


m2 


w2 


2009-July-25 


55037 


> 18.78 


19.31±0.28 


18.68±0.23 


18.57±0.18 


18.25±0.16 


18.49±0.12 


2012-Apr-25 


56042 


18.65±0.33 


> 19.80 


18.98±0.27 


18.98±0.22 


18.85±0.10 


18.58±0.11 


2012-May-25 


56072 


18.72±0.31 


19.28±0.22 


18.74±0.06 


18.32±0.14 


18.60±0.17 


18.35±0.10 


2012-June-25 


56103 


> 18.76 


19.12±0.27 


18.67±0.27 


18.51±0.08 


18.25±0.15 


18.62±0.13 


2012-Aug-07/08 


56146/7 


19.39±0.30 


20.12±0.28 


18.90±0.14 


18.79±0.10 


18.61±0.08 


18.65±0.06 
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Figure 2. 15 GHz radio light curve of PKS 1502+036 for the pe- 
riod 2008 August 7-2012 November 4 from the OVRO telescope. 



care; the solution interval (30 min) was chosen to be longer 
than the scan-length to remove residual systematic errors 
and fine tune the flux density scale, but not to force the in- 
dividual data points to follow the model. 
Additional archival data have been considered to improve 
the time sampling. A priori amplitude calibration was de- 
rived using measurements of the system temperature and 
the antenna gains. Uncertainties of the amplitude calibra- 
tions were assumed to be 5%. For MDJAVE data, we i mported 
the fully calibrated uv data sets (|Lister et al.ll2009l ). 
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Figure 3. Multi-epoch radio spectra of PKS 1502-1-036 observed 
with VLA and VLBA during the monitoring campaign of the 
'bright' HFP sample. Epoch a: X symbol and solid line; epoch h: 
empty circle and dot line; epoch c: empty star and dot-dashed 
line; epoch I (VLBA): filled circle and dash-and-three-dot line; 
epoch d: triangle and dashed line. 



flux density scale is calibrated using obs ervations of 3 C 286 , 
adopting the flux density (3.44 Jy) from lBaars et al] |l973). 
This results in about a 5% absolute scale uncertainty, which 
is not reflected in the plotted errors. As shown in Fig. [S] 
PKS 1502-1-036 was highly variable at 15 GHz during the 
OVRO monitoring, with a flux density spanning 376 mjy 
(at MJD 55726) to 669 mJy (at MJD 56139). 



5.2 OVRO 

As part of an ongoing blazar monitoring programme, the 
OVRO 40 m radio telescope has observed PKS 1502-1-036 
at 15 GHz regularly since the end of 2007 l|Richards et al.l 
l201lh . This monitoring programme includes over 1500 
known and likely 7-ray loud blazars above declination —20°. 
The sources in this programme are observed in total in- 
tensity twice per week with a 4 mJy (minimum) and 3% 
(typical) uncertainty in their flux densities. Observations 
are performed with a dual-beam (each 2.5 arcmin FWHM) 
Dicke-switched system using cold sky in the off-source beam 
as the reference. Additionally, the source is switched be- 
tween beams to reduce atmospheric variations. The absolute 



5.3 Radio data: results 

5.3.1 The radio spectrum 

Simultaneous multifrequency VLA observations carried out 
at various epochs showed substantial spectral and flux den- 
sity variability. To determine th e change in the spe ctral 
peak, we follow the approach by lOrienti et al.l (|2010l ) and 
we fit the simultaneous radio spectrum at each epoch with 
a purely analytical function: 

log S = a -I- log X (b -f c log v) (1) 

where S is the flux density, v the frequency, and a, b and c 
are numeric parameters without any direct physical mean- 
ing. The best flts to the spectra are shown in Fig. |3l and the 
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Table 3. Logs of VLA and multifrequency VLBA observations. Column 1: observation date in UT; Column 2: observation date in MJD; 
C olumn 3: observation c ode; Column 4: array; Col umn 5: observinR b ands; Column 6: c ode used in this pap er; Column 7: reference: 
l= lDallacasa et al.l ll200d) ; 2= lTinti et al] ^M); 3= IOrienti et alj ||2007| '); 4=this paper; 5= IOrienti et al.l l l2oi36t) . 



Date (UT) 


MJD 


Observation code 


Array 


Bands 


Code 


Reference 


1999-Sep-25 


51446 


AD428 


VLA-A 


L,C',X,U,K 


a 


1 


2002-July-03 


52458 


AT275 


VLA-B 


L,C,X,U,K 


b 


2 


2003-Sep-13 


52895 


AD488 


VLA-AnB 


L,C,X,U,K 


c 


3 


2007-Apr-14 


54204 


AO210 


VLA-D 


L,C,X,K 


d 


4 


2001-Oct-06 


52188 


BF068 


VLBA'' 


U 


e 


4 


2002-Jan-ll 


52285 


BD077 


VLBA 


X,U 


f 


5 


2002-Sep-26 


52481 


RDV35 


VLBA 


X 


9 


4 


2005-Dec-07 


53711 


BF075 


VLBA'' 


u 


h 


4 


2006-June-26 


53912 


BF075 


VLBA'' 


u 


i 


4 


2006-July-22 


53938 


BO025 


VLBA 


L,S,C,X,U 


I 


4 


2010-Oct-15 


55484 


BL149 


VLBA^ 


U 


m 


4 


2010-Oct-25 


55494 


BL149 


VLBA^ 


U 


n 


4 


2011-May-06 


55687 


BL149 


VLBA'' 


U 





4 


2012-Jan-02 


55928 


BL178 


VLBA^ 


U 


P 


4 


2012-Sep-27 


56197 


BL178 


VLBA'' 


u 


g 


4 



"archival VLBA data; ''data from the MQJAVE programme 



Table 4. VLA flux density. Column 1: code from Table |3] Columns 2-9: flux density at 1.4, 1.7, 4.5, 4.9, 8.1, 8.4, 15.3, 22.2 GHz, 
respectively; Column 10: peak frequency; Column 11: peak flux density; Columns 12 and 13: spectral index below and above the peak 
frequency, respectively; Column 14: variability index (see Section l5.3.ip . 



Code 


■51.4 




54.5 


S'4.9 


58.1 


58.4 


5l5.3 


522.2 




5p 




Ola 


V 




(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(GHz) 


(mJy) 








a 


455 


547 


921 


929 


901 


895 


830 


738 


7.7±0.8 


940 


-1.7 


0.2 


107 


b 


413 


465 


740 


741 


721 


710 


665 


567 


7.2±0.7 


745 


-0.5 


0.2 


3 


c 


382 


428 


608 


620 


620 


610 


515 


468 


6.5±0.6 


622 


-0.4 


0.3 


122 


d 


332 




679 


692 


720 


720 




611 


8.8±0.9 


737 


-0.5 


0.2 


38 



derived peak frequencies at the various epochs are reported 
in Tables |4] and [5] Statistical errors derived from the fit are 
not representative of the real uncertainty on the estimate 
of the peak frequency. For this reason, we prefer to assume 
a conservative uncertainty on the peak frequency of 10%. 
In addition, we computed the spectral index in the optically 
thick (ob) and thin part (oa) of the spectrum. Both the spec- 
tral shape and the spectral peak change randomly among the 
various observing epochs, similar to what is found in blazars 
(e.g. _Torniairien ct al. 2007.) as well as in radio-lo ud NLSyls 
(e.g. lYuan et al.li2008l : iD'Ammando et al.ll2012al). 

Following the approach of lOrienti et al. I l|2007l ). we es- 
timated the spectral variability by means of the variability 
index: 



V = 



in 



8 = 1 



(2) 



where Si is the flux density at the ith frequency mea- 
sured at one epoch. Si is the mean value of the flux density 
computed by averaging the flux density at the ith frequency 
measured at all the available epochs, Oi is the root mean 
square (rms) on Si — Si and m is the number of sampled 
frequencies. The variability index computed at the differ- 
ent epochs ranges between 3 and 122, indicating that this 
source is highly variable (Table [l]) . We prefer to compute 



the variability index for each epoch instead of considering 
all the epochs together in order to potentially detect small 
outbursts. 



5.3.2 Radio structure 

PKS 1502-1-036 is unresolved on arcsecond scales typical of 
the VLA (Fig. 0]). However, when imaged with the parsec 
scale resolution provided by VLBA observations, its radio 
structure is marginally resolved and a second component 
seems to emerge from the core. In the VLBA images at 15 
GHz, the radio structure is clearly extended, suggesting a 
core-jet-like morphology (Figs [S] [6] and [7]). The radio emis- 
sion is dominated by the core component, which is unre- 
solved in the 15 GHz image with an angular size < 0.4 mas. 
The jet-like feature accounts for about 4% of the total flux 
density in the observation presented here, and it extends to 
about 3 mas (Fig. [7]). Since the VLBA flux density is in good 
agreement with the nearly contemporaneous measurements 
by the OVRO 40 m single dish, we can exclude significant 
emission from extended, low surface brightness features. 

To investigate possible changes in the radio structure, 
we modelled the performed visibility-based model fitting us- 
ing the DIFMAP software. The radio structure consists of a 
bright compact core plus weak extended emissions westward 
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Table 5. Multifrequency VLBA flux density. Column 1: code from Tabled Columns 2-6: flux density at 1.7, 2.3, 5.0, 8.4, 15.3 GHz, 
respectively; Columns 7 and 8: peak frequency and peak flux density, respectively (see Section l5.3.1|l : Columns 9 and 10: spectral index 
below and above the peak frequency, respectively. 



Code 51.7 


■52.3 


■5*5.0 


■58.4 


■5l5.3 




Sp 




(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(GHz) 


(mJy) 




/ 






630 


566 








9 






800 










I 365 


480 


670 


701 


547 


6.6±0.7 


704 


-0.6 0.4 



and is well fitted with a point source (core) plus two (cir- 
cular) Gaussian components ('CI' and 'C2'); we could not 
identify the third component for the data of 2006 March 20 
due to limited dynamic range. The results of image anal- 
yses are listed in Table |B] The error in total flux density 
was determined from the uncertainties of the amplitude cal- 
ibration and image rms noise. The error in the measured 
separation between the core and the jet components was de- 
termined from the root-sum-square of the uncertainties of 
the position measurements for CI (or C2) and the core. The 
uncertainty in the relative positions is estimated from 

Ar = a/(S'p/rms), (3) 

where a is the major axis of the deconvolved component, 
Sp is its peak flux density an d rms is the la noise leve l 
measured on the image plane ijPolatidis fc Conwavlliooi ). 
In case the errors estimated by equation ((3| are unreliably 
small, we assume a more conservative value for Ar that is 
10% of the FWHM of the beam. 

In contrast to what is found in the radio-loud NLSyl 
SBS 0846-1-513, where an ap parent superluminal motion has 
been detected {P ~ 10.9c. iD'Ammando et"al] 12012^ 1 . no 
signiflcant proper motion is detected for the jet components 
of PKS 1502+036 (Fig.©. 



PKS 1502+036 22.2 GHz VLA 




15 10 5 -5 -10 -15 

ARC SEC 



Figure 4. VLA image at 22.2 GHz of PKS 1502-f036. On the 
image, we provide the restoring beam, plotted in the bottom-left 
corner, the peak flux density in mjy/beam and the flrst contour 
(f.c.) intensity in mJy/beam, which is three times the off-source 
noise level. Contour levels increase by a factor of 2. 



6 DISCUSSION AND CONCLUSIONS 

The detection of four radio-loud NLSyl s in 7-rays during 
the first year of Fermi-LAT operation jAbdo et all 120091 ) 
has stimulated the interest in this class of AGN. Currently 
the search for new 7-ray emitting NLSyl s has led to only one 
new so urce detected at high significance (|D'Ammando et al.l 
l2012al ). So far most of the attention has been directed to- 
wards the two NLSyls that showed 7 -ray flaring activity , 
PMN J0948-H0022 and SB S 0846-^513 l|Foschini et al.ll201ll : 
ID'Ammando et al.|[2012ah . 

PKS 1502-1-036 showed no signiflcant 7-ray variability 
during 2008 August 4-2012 November 4, with the 0.1-100 
GeV flux ranging between (3-7)xl0~* photons cm~^ s^^ 
using 3-month time bins. The average LAT spectrum accu- 
mulated over 51 months is well described by a power-law 
with a photon index of F = 2.60 ± 0.06, which is similar 
to the mean value observed for flat spectrum radio quasars 
(FSRQs) during t he flrst two years of Ferm i-hNT operation 
(F = 2.42 ± 0.17: lAckermann et al.ll2012al '). 
In the same way, the average apparent isotropic 7-ray lu- 
minosity of PKS 1502+036 is 7.d,x\(f^ erg s'^ in the 
0.1-100 GeV range, a typical value for a FSRQ (e.g. 



iGrandi fc Torresill2012l ). and it is in agreement with what 
was found fo r SBS 0846+513 and PM N J0948+0022 (see 
discussion in ID'Ammando et al.|[2012al ). This suggests the 
presence of a relativistic jet with Doppler factors (S) as 
large as in blazars. Modelling of the spectral energy distri- 
bution (SED) of the radio-loud NLSyls detecte d by Fermi- 
LAT indicates Doppler factors larger than 10 (|Abdo et al.l 
I2OO9I : ID'Ammando et al.llioiial 'l. and in particular 5 ~ 18 
for PKS 1502+036. However, the study of the proper mo- 
tion could not detect any signiflcant motion for this source, 
while the radio spectral variability and the one-sided struc- 
ture further require the presence of boosting effects in a 
relativistic jet. 

The 15 GHz light curve (Sec. I5.2|l is highly variable. 
In particular, two outburst episodes are clearly present: the 
flrst peaks at MJD 54923, and the second peaks at MJD 
56139. These outbursts seem to take place when a flux en- 
hancement is also observed in 7-rays. Unfortunately, the 
statistics in 7-rays are not adequate to allow a study of the 
variability on time-scales short enough to be related to the 
radio monitoring, while in the optical and UV bands the 
time sampling was too poor. For this reason, we estimate 
the variability time-scale At on the basis of the radio data, 
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Table 6. Log and results of VLBA observations at 15 GHz. 



Date (UT) MJD Restoring beam Image noise VLBA total flux CI position C2 position 
(mas X mas) (mjy/beam) (mjy) (mas) (mas) 



2001-Oct-06 


52188 


1.01 X 0.41 


2.2 


439 ± 22 


0.29 ± 0.07 


1.04 ±0.17 


2002-Jan-ll 


52285 


1.06 X 0.44 


5.2 


565 ± 29 


0.23 ± 0.07 


1.21 ±0.56 


2006-Mar-20 


53814 


1.58 X 0.62 


2.0 


608 ± 30 


0.38 ± 0.10 




2006-June-27 


53913 


1.08 X 0.36 


1.1 


553 ± 28 


0.18 ± 0.06 


1.08 ±0.13 


2006-July-21 


53937 


1.60 X 0.71 


0.8 


547 ± 27 


0.22 ± 0.11 


1.13 ±0.11 


2010-Oct-15 


55484 


1.22 X 0.51 


0.2 


504 ± 25 


0.30 ±0.08 


1.00 ±0.09 


2010-Oct-25 


55494 


1.33 X 0.52 


0.2 


487 ± 24 


0.31 ± 0.08 


0.97 ±0.10 


2011-May-26 


55707 


1.32 X 0.63 


0.3 


429 ± 21 


0.32 ± 0.09 


0.98 ±0.13 


2012-Jan-02 


55928 


1.26 X 0.53 


0.2 


444 ± 22 


0.33 ± 0.08 


0.94 ±0.09 


2012-Sep-27 


56197 


1.23 X 0.53 


0.3 


588 ± 29 


0.17 ±0.08 


0.90 ±0.12 




Figure 5. VLBA image at 15.3 GHz of PKS 1502+036 collected 
on 2002 January 11. On the image, we provide the restoring 
beam, plotted in the bottom-left corner, the peak flux density 
in mJy/beam and the first contour (f.c.) intensity in mJy/beam, 
which is three times the off-source noise level. Contour levels in- 
crease by a factor of 2. 



Figure 6. VLBA image at 15.3 GHz of PKS 1502+036 col- 
lected on 2006 July 21. On the image, we provide the restoring 
beam, plotted in the bottom-left corner, the peak flux density 
in mJy/beam and the first contour (f.c.) intensity in mJy/beam, 
which is three times the off-source noise level. Contour levels in- 
crease by a factor of 2. 



similarly to IValtoia et"aD (|l999ll . For At, we assume the 
time interval of the flux density variation to be between the 
minimum and maximum flux density of a single outburst 
lAS'l. This assumption implies that the minimum flux den- 
sity corresponds to a stationary underlying component and 
the variation is due to a transient component. Taking into 
consideration the time dilation due to the cosmological red- 
shift we flnd that the intrinsic time lag is Ar = At/(1 + z), 
while the intrinsic flux density variation at the observed fre- 
quency is |ASi| = |AS| X (1 + 2)^"°. The brightness tem- 
perature at a given frequency is 



1 S{u) 
2k Q 



(4) 



where k is the Boltzmann constant, Q, is the solid angle of 
the emitting region, and c is the speed of light. Using the 



causality principle, we can determine the angular size of the 
region responsible for the outburst 

cAt {1 + zf 
(l + z) Dr. ■ 
Knowing that the solid angle is 

we derive the rest-frame brightness temperature from 

, ^ 2 \AS\Dl 
^ nk Af2iy'2{l + zy+" ■ 

During the two outburst episodes we have flux density vari- 
ability AS — 185 and 216 mJy, respectively, while At = 75 
and 85 d, respectively. If in equation (0 we consider these 
values, and we assume a = 0.3, i.e. the average value ob- 
tained by fitting the optically thin spectrum (Tables |4] and 



(6) 



(7) 
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Figure 7. VLBA image at 15.3 GHz of PKS 1502+036 collected 
on 2010 October 15. On the image, we provide the restoring 
beam, plotted in the bottom-left corner, the peak flux density 
in mJy/beam and the first contour (f.c.) intensity in mJy/beam, 
which is three times the off-source noise level. Contour levels in- 
crease by a factor of 2. 



2002 



2004 



2006 2008 2010 2012 




Figure 8. Angular separation between the jet components (CI 
and C2, full and empty circles respectively) and the core in 
PKS 1502-1-036. 



[5}, we obtain T'g ~ 2.5x10^"^ K, which exceeds the value 
derived for the Compton catastrophe. Assuming that such a 
high value is due to Doppler boosting, we can estimate the 
variability Doppler factor 5var, by means of 

l/(3+a) 

S.ar = ( 7^ I , (8) 



Tint 

where Tint is the intrinsic brightness temperature. Assuming 



a typi cal value Ti„t = 5x 1 K , as derived by e.g.lReadhead 
' 19941) : iHovatta et"al] (|2009l '): iLahteenmaki fc Vahaoia 
19991 ). we obtain S^ar = 6.6. 



Another way to derive a lower limit to the Doppler 
factor is by means of the jet/counter-jet brightnesses ra- 
tio. Assuming that both the jet and counter-jet have the 
same intrinsic power, the different brightness observed can 
be related to Doppler boosting by means of 



1 + ficose \ ' 
/3cos6l ) 



1 



(9) 



where Bj and Bcj are the jet and counter-jet brightness, 
respectively, Q is the viewing angle, /3 is the bulk velocity in 
terms of the speed of light. We prefer to compare the surface 
brightness instead of the flux density because the jet has a 
smooth structure without clear knots. We measured the jet 
brightness as 20 mJy/beam in the 15 GHz image. In the 
case of the counter-jet, which is not visible, we assumed an 
upper limit for the surface brightness that corresponds to 
0.15 mJy/beam, i.e. 1(t noise level measured on the image. 
From the brightness ratio estimated from equation ((9l) we 
obtain /3cos& > 0.8, implying that the minimum velocity is 
/3 > 0.8 (corresponding to a minimum bulk Lorentz factor 
r = 2.8) and a maximum viewing angle Q = 36°. From the 
minimum values derived, we can estimate a lower limit to 
the Doppler factor by means of 

' - r(i-'/3cos^) (1°) 

and we obtain S > 1.7. We remark that the Doppler factors 
estimated by means of variability and jet/counter-jet bright- 
ness ratio a re lower limits. The Doppler factor obtained in 
lAbdo et ai] ([2009) by modelling the SED is much larger, S 
= 18. The discrepancy may be a consequence of the different 
emitting regions. The radio variability is likely produced in 
the core component instead of along the jet. Furthermore, 
the 7-ray emission can be produced in a very compact re- 
gion, as suggested by the ver y short variability tim e-scale 
observed in some FSRQs (e.g. iTavecchio et al.ll201ol ). Such 
a small region is usually self-absorbed at the typical radio 
frequencies (cm and mm wavelen g ths). Indeed, the fit used 
to model the SED in lAbdo et al.1 (|2009f ) does not take into 
consideration the radio data. On the other hand, the region 
responsible for the radio emission is related to regions of the 
jets that are further away from the central AGN, where the 
opacity is less severe. 

A black hole mass of 4x10*^ Mq was estimated by 
lYuan et al.l (|2008l ) for PKS 1502-1-0 36 on the ba.sis of the 
H^g broad line. On the other hand, lAbdo et all (|2009l ) de- 
rived a black hole mass of 2x10^^ Mq by means of the SED 
modelling of the source. As in the case of SBS 0846+513, ra- 
dio and 7-ray properties of PKS 1502+036 seem to show the 
characteristics of a (possibly young) blazar at the low end 
of the black hole mass distribution. A core-jet structure has 
been resolved at 15 GHz in VLBA images of PKS 1502+036. 
The main difference with respect to SBS 0846+513 seems to 
be a fainter jet-like structure and no observed apparent su- 
perluminal motion. 

In Fig. |9] we compare the 7-ray light curve collected 
by Ferrm-LAT during 2008 August-2012 November to the 
X-ray (0.3-10 keV), UV (m2 fiher), optical {u and V fihers) 
and radio (15 GHz) light curves collected by Swift, CRTS 
and OVRO. No significant flaring activity was detected in 
7-rays. The lack of significant 7-ray variability is not sur- 
prising taking into account that only 40% of the AGN in 
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Figure 9. Multifrcqucncy light curve for PKS 1502+036. The period covered is 2008 August-2012 November. The data sets were 



collected (from top to bottom) by Fermi-LAT (7-rays), Swi/f-XRT (0.3-10 keV, in units of 10 ^'^ org cm 
u bands, in units of mjy), CRTS (V band, in units of mjy) and OVRO (15 GHz, in units of mjy). 



-1), Swift-VYOT {w2 and 



the second LAT AGN C ataloRue (2LAC) clean s ample has 
shown a clear variability l|Ackermann et al]|2012al ). The con- 
tinuous monitoring of the 7-ray sky provided by f ermi-LAT 
will allow us to catch, if it happens, a flaring activity from 
PKS 1502+036. Interestingly a slight increase from radio to 
UV was observed at the end of 2012 June during a period 
of relatively high 7-ray flux. Similarly, the maximum flux 
density in V band was observed when the 7-ray emission 
was increasing. However, the sparse coverage in optical-to- 
X-rays does not allow us to obtain conclusive evidence. A 
regular monitoring of this source, as well as the other 7-ray 
emitting NLSyls, from radio to 7-rays will be fundamental 
for continuing to investigate the nature and the emission 
mechanisms of these objects. This will be crucial for reveal- 
ing differences and similarities between 7-ray NLSyls and 
blazars. 
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